Finding the unloaded geometryΒΆ
from fenics_plotly import plot
import pulse
geometry = pulse.HeartGeometry.from_file(pulse.mesh_paths["simple_ellipsoid"])
# geometry = pulse.geometries.prolate_ellipsoid_geometry(mesh_size_factor=3.0)
material = pulse.NeoHookean()
# material = pulse.Guccione()
2022-01-21 12:55:09,141 [1095] INFO pulse.geometry_utils:
Load mesh from h5
# Parameter for the cardiac boundary conditions
bcs_parameters = pulse.MechanicsProblem.default_bcs_parameters()
bcs_parameters["base_spring"] = 1.0
bcs_parameters["base_bc"] = "fix_x"
# Create the problem
problem = pulse.MechanicsProblem(geometry, material, bcs_parameters=bcs_parameters)
# Suppose geometry is loaded with a pressure of 1 kPa
# and create the unloader
unloader = pulse.FixedPointUnloader(problem=problem, pressure=3.0)
2022-01-21 12:55:09,350 [1095] INFO pulse.unloader:
-------------------------- Unloading options ---------------------------
Target pressure: 3.0
maxiter = 10
tolerance = 0.0001
regenerate_fibers (serial only)= False
------------------------------------------------------------------------
# Unload the geometry
unloader.unload()
2022-01-21 12:55:09,361 [1095] INFO pulse.unloader: ------------------------------------------------------------------------
2022-01-21 12:55:09,362 [1095] INFO pulse.unloader: ----------------------------Start unloading-----------------------------
2022-01-21 12:55:09,364 [1095] INFO pulse.unloader: ------------------------------------------------------------------------
2022-01-21 12:55:09,369 [1095] INFO pulse.unloader:
LV Volume of original geometry = 0.749 ml
2022-01-21 12:55:09,477 [1095] INFO pulse.unloader:
Do an intial solve
2022-01-21 12:56:38,316 [1095] INFO pulse.unloader:
Iteration: 0
Ignoring precision in integral metadata compiled using quadrature representation. Not implemented.
2022-01-21 12:56:44,393 [1095] WARNING FFC: Ignoring precision in integral metadata compiled using quadrature representation. Not implemented.
Ignoring precision in integral metadata compiled using quadrature representation. Not implemented.
2022-01-21 12:56:51,514 [1095] WARNING FFC: Ignoring precision in integral metadata compiled using quadrature representation. Not implemented.
2022-01-21 12:56:52,384 [1095] INFO pulse.unloading_utils:
LV Volume of original geometry = 0.568 ml
2022-01-21 12:58:21,653 [1095] INFO pulse.unloading_utils:
LV Volume of inflated geometry = 0.568 ml
2022-01-21 12:58:22,134 [1095] INFO pulse.unloading_utils:
LV Volume of new reference geometry = 0.717 ml
2022-01-21 12:58:22,177 [1095] INFO pulse.unloader:
Residual: 0.011328623428538229
2022-01-21 12:58:22,178 [1095] INFO pulse.unloader:
Iteration: 1
2022-01-21 12:58:22,653 [1095] INFO pulse.unloading_utils:
LV Volume of original geometry = 0.596 ml
2022-01-21 12:58:27,860 [1095] INFO pulse.unloading_utils:
LV Volume of inflated geometry = 0.596 ml
2022-01-21 12:58:28,344 [1095] INFO pulse.unloading_utils:
LV Volume of new reference geometry = 0.754 ml
2022-01-21 12:58:28,390 [1095] INFO pulse.unloader:
Residual: 0.0016980102145402018
2022-01-21 12:58:28,393 [1095] INFO pulse.unloader:
Iteration: 2
2022-01-21 12:58:28,912 [1095] INFO pulse.unloading_utils:
LV Volume of original geometry = 0.592 ml
2022-01-21 12:58:34,027 [1095] INFO pulse.unloading_utils:
LV Volume of inflated geometry = 0.592 ml
2022-01-21 12:58:34,512 [1095] INFO pulse.unloading_utils:
LV Volume of new reference geometry = 0.749 ml
2022-01-21 12:58:34,558 [1095] INFO pulse.unloader:
Residual: 0.0002592876931297112
2022-01-21 12:58:34,561 [1095] INFO pulse.unloader:
Iteration: 3
2022-01-21 12:58:35,052 [1095] INFO pulse.unloading_utils:
LV Volume of original geometry = 0.593 ml
2022-01-21 12:58:40,194 [1095] INFO pulse.unloading_utils:
LV Volume of inflated geometry = 0.593 ml
2022-01-21 12:58:40,680 [1095] INFO pulse.unloading_utils:
LV Volume of new reference geometry = 0.750 ml
2022-01-21 12:58:40,723 [1095] INFO pulse.unloader:
Residual: 3.9313042959387195e-05
2022-01-21 12:58:40,726 [1095] INFO pulse.unloader: ########################################################################
Unloading suceeding
# Get the unloaded geometry
unloaded_geometry = unloader.unloaded_geometry
fig = plot(geometry.mesh, opacity=0.0, show=False)
fig.add_plot(plot(unloaded_geometry.mesh, color="red", show=False))
fig.show()